Ultracold atoms in an optical lattice with dynamically variable periodicity 
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The use of a dynamic "accordion" lattice with ultracold atoms is demonstrated. Ultracold atoms 
of ^^Rb are trapped in a two-dimensional optical lattice, and the spacing of the lattice is then 
increased in both directions from 2.2 /xm to 5.5 /im. Atoms remain bound for expansion times as 
short as a few milliseconds, and the experimentally measured minimum ramp time is found to agree 
well with numerical calculations. This technique allows an experiment such as quantum simulations 
to be performed with a lattice spacing smaller than the resolution limit of the imaging system, while 
allowing imaging of the atoms at individual lattice sites by subsequent expansion of the optical 
lattice. 



PACS numbers: 37.10.Jk, 03.75.Lm, 03.67.-a 

Optical lattices create clean, tunable and flexible peri- 
odic potentials for ultracold atoms that are an important 
tool for investigating the quantum behavior of strongly 
correlated many-body systems [1^3] . Such investigations 
require lattices with sub-micron spacing to ensure that 
the quantum dynamics occurs on millisecond timescales 
(and so is not subject to decoherence). The small spac- 
ing also ensures that the tunneling energy J and on-site 
interaction energy U are comparable ( J//7 ~ 1) at a lat- 
tice depth for which the band gap Eg ^ /c^T, assuming 
a typical temperature of tens of nanokelvin. The sub- 
micron lattice spacing makes it challenging to observe 
atoms at individual lattice sites with light or near in- 
frared radiation, and typically time- of- flight expansion 
has been used to probe the momentum of the atoms [4| . 
Being able to detect the position of the atoms in situ 
can yield further information about the system, and the 
ability to address atoms at single lattice sites is crucial 
for quantum information processing [5]. In situ imaging 
of the atoms has been demonstrated for lattice spacings 
from 2 jivd upwards [6. .8,, however at these spacings the 
single-atom tunneling is negligible on typical experimen- 
tal timescales. Scanning electron microscopy has been 
used to detect single atoms in a 0.6 jivn optical lattice 
and determine the density distribution by summing over 
multiple images [9^. Recently there has been significant 
progress on directly imaging single atoms trapped in op- 
tical lattice sites at sub-micron lattice scales using flu- 
orescence imaging [TOl [TJ: in both cases sophisticated 
optical arrangements were used. The difficulty of ob- 
taining the theoretical maximum resolution in an optical 
imaging system increases dramatically as the numerical 
aperture increases. 

Our method of imaging such systems is to use a lattice 
with a dynamically variable spacing. Ultracold atoms 
are prepared in a lattice with small spacing and then the 
lattice potential is expanded to facilitate imaging. This 
allows individual lattice sites to be resolved without the 
need for a complex or expensive microscope arrangement. 
This approach could also be used for addressing individ- 
ual lattice sites; a previous experiment [8 demonstrated 



addressing single lattice sites spaced by 5.3 fivn. Increas- 
ing the spacing between the trapped atoms also reduces 
the error in qubit readout fidelity introduced by cross 
talk between neighboring atoms [T2^. Dynamically vary- 
ing the lattice spacing could provide a promising tech- 
nique for measurement based quantum computing where 
very high ground state fidelities are required — after a 
cluster state is created [13l [14] the optical lattice could 
be expanded to allow measurements to be carried out. 

Diffraction from a ID lattice with variable periodicity 
[15] , and investigation of the dynamics of a Bose-Einstein 
condensate (BEG) loaded into a large period ID lattice 
with tuneable spacing [16 have been demonstrated; nei- 
ther of these arrangements, however, were capable of 
changing the lattice spacing whilst keeping the atoms 
trapped. An arrangement for generating ID dynamically 
varying accordion lattices has been shown [17] in which 
mechanical means were used to vary the lattice spacing 
from 0.96 jiia to 11.2 jiia over one second (although it 
was not demonstrated with trapped ultracold atoms). 

Here we demonstrate the first use of an accordion lat- 
tice arrangement to dynamically increase the optical lat- 
tice spacing while the atoms are trapped. Unlike in the 
work described above [ iSHTT] , the lattice spacing is in- 
creased in two dimensions, and the arrangement allows 
the optical potential to be varied smoothly and quickly, 
free of mechanical vibrations. We show that the atoms 
remain trapped in the optical lattice for ramp times as 
short as a few milliseconds, well within the timescales of 
a typical BEG experiment. 

The accordion lattice is generated using dual axis 
acousto-optic deflectors (AODs), and an optical arrange- 
ment described in detail in [18 . Briefly, the optical lat- 
tice is formed in the focal plane of a multi-element objec- 
tive lens made according to the design described in [19]. 
The intersection angle of the beams determines the lat- 
tice spacing d, i.e. d = XF/D^ where A is the wavelength 
and F is the focal length of the lens. D is the distance 
between the beams in the back focal plane of the lens, 
which is controlled by the deflection angle from an AOD. 
Thus the lattice spacing can be varied dynamically by 



changing the frequency of the radio frequency radiation 
controhing an AOD. This arrangement provides an ex- 
tremely flexible optical lattice which can be rotated [20] , 
or expanded and contracted as described here. The mini- 
mum lattice spacing achievable is given by the numerical 
aperture of the lens, which in our case is 0.27. The opti- 
cal lattice is formed by the interference of near infrared 
radiation at 830 nm provided by a Ti:sapphire laser (Co- 
herent MBR-110). The optical lattice potential has the 
form 



V{x,y,t) = Vi{x,y){cos^[nx/d{t)] 



cos^ 



[7ry/d{t)]}; (1) 



the spatially-dependent lattice depth Vi{x, y) has a Gaus- 
sian distribution — Vbexp[— 2(x^ + y'^)/w'^], where Vq is 
the peak lattice depth and w is the waist. 

We start with a nearly pure BEC of 80,000 ^^Rb atoms 
in the F = l^mp = —1 state in a combined mag- 
netic and optical trap with harmonic trapping frequen- 
cies {ujr^^z} = 27r X {20.1, 112} Hz. This trap is on for 
the remainder of the experimental sequence. An opti- 
cal lattice potential of the form given in Eq. (fTl) with 
Vq = h X 1.5 kHz, w = 69 /im and d = 2.2 /im is ramped 
on in 500 ms. This timescale is long enough to ensure 
the system is in the ground state of the optical lattice. 
Figure 1(a) shows an in situ absorption image of atoms 
trapped in the optical lattice; we first optically pump the 
atoms to the F = 2 hyperfine level, and then probe the 
atoms by applying a 20 /is pulse of circularly polarized 
light resonant with the 56'i/2 F = 2 - 5P3/2 F^ = 3 tran- 
sition. The optical pumping light is detuned to ensure 
only a fraction of the atoms are pumped to the F = 2 
level, reducing the optical density in the lattice to order 
one to prevent saturation of the absorption images. The 
individual lattice sites cannot be resolved, because the fi- 
nite axial extent of the cloud, which is 11 /im diameter or 
more, is comparable with the depth of field of the imag- 
ing system. The trapped cloud has radius ro = 18 /im, 
which indicates that there are around 8 lattice sites filled 
in each direction, with around 1200 atoms in each central 
lattice site. 

The lattice depth is then ramped to /i x 37 kHz in 
5 ms, which cuts off tunneling between the lattice sites. 
This ensures that the lattice dynamics are frozen during 
the expansion of the lattice. The ramp time is short 
compared to the single-particle tunneling time, however 
it is long enough to ensure that the atoms remain in the 
lowest energy band of the optical lattice, i.e. {dH/dt) <C 

The lattice spacing is varied as a smoothed ramp 
given by the Gauss error function d{t) = do{l -\- {S — 
l){eTf[7.U{t/tr - 0.5)] + l}/2), where do is the initial 
spacing (2.2 /im), S is the ratio of the final spacing to 
the initial spacing and t^ is the ramp time (which is var- 
ied). The peak lattice depth Vq is kept constant during 
the ramp, and 5* is 2.5. In principle the lattice spacing 
can be expanded to any size provided that the lattice 
spacing remains smaller than the beam waist w. In prac- 
tice, to ensure that atoms at the edge of the cloud are not 



lost, S < w/ro should be satisfied. For our parameters, 
after a 2.5 times expansion, the radius is 45 /im, and the 
depth at the outermost lattice site has dropped to about 
40% of the depth at the center of the lattice. The final 
lattice spacing, 5.5 /im, is large enough to be resolved by 
our imaging system, as shown in Fig. |l(b) The num- 
ber of lattice sites filled is consistent with the dimensions 
of the cloud before expansion. To confirm that the ac- 
cordion lattice is working as expected, we also expand 
by a factor of 2.5 from an initial lattice spacing that is 
large enough to be resolved before expans ion ( 5.5 / im to 
13.75 /im). The result, presented in Fig. |l(c)] and p^d)| 
shows that the atoms do remain bound for this ramp. 
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(c) Before ramp 



(b) After ramp 
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(d) After ramp 



FIG. 1. Absorption image of atoms trapped in a two- 
dimensional optical lattice before and after the lattice spacing 
has been ramped by: a factor of 2.5 from 2.2 /xm to 5.5 /im in 
10 ms [(a) and (b)]; a factor of 2.5 from 5.5 /xm to 13.75 /im 
in 50 ms [(c) and (d)]. The field of view is 120 /xm by 120 /xm 
in each image. 



During the expansion of the lattice, the central lattice 
site remains stationary [21 , while the first lattice site 
moves with speed yi = d, the second lattice site moves 
with speed V2 = 2d, and the n*^ lattice site moves with 
speed Vn = nd. To determine the evolution of the system 
during the ramp in lattice spacing, we describe the atoms 
at each lattice site by localized wave functions IV^nc^n^), 
evolving according to 



ih 



dt 



\2M 



+ T/(x,y,i) |^„^. (i)), (2) 



where M is the mass of an atom 
ignored. 

If the accordion lattice is to be used for quantum com- 
puting, it is important that the ground state fidehty re- 
mains high during the lattice expansion. However, if the 
lattice is expanded for the purpose of imaging the atoms, 
moving the atoms slowly enough to retain high ground 
state fidelity is not necessary — instead there is the less 
stringent requirement that an atom remains in a bound 
state of the same lattice site. 

To estimate the maximum acceleration for which the 
atoms in the n*^ well remain bound, we can move to 
an accelerating reference frame in which the n*^ lattice 
site is at rest. In this reference frame, the lattice poten- 
tial is tilted due to a linear potential term Mndx [22] 
which has the effect of reducing the potential barrier be- 
tween neighboring lattice sites, as illustrated in Fig. |2J 
We would expect atoms to be lost when the minimum 
barrier height /SV approaches the energy of the trapped 
atoms, which for our experimental parameters would oc- 
cur for accelerations corresponding to a ramp time of 
around 3 ms. /SV decreases with increasing scaled accel- 
eration A = Mndd/TrVo^ and so the maximum accelera- 
tion can be larger for a deeper lattice, or for a smaller 
lattice spacing. This approach gives an order of magni- 
tude estimate for the minimum expansion time, however 
it does not take into account the reduction of the barrier 
height during expansion due to the envelope of the lat- 
tice beams, the effect of the magnetic trap, nor does it 
take into account excitations to higher bands during the 
ramp. 



and interactions are along the direction of motion for each lattice site 
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FIG. 2. Effective trapping potential at the n* lattice site in 
a reference frame accelerating at nd, in terms of the scaled 
acceleration A = Mndd/TrVo. 
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FIG. 3. (color online) (a) Numerically calculated probabil- 
ity of finding a single particle in the ground state of the 
cos^(x) potential well of the n^ lattice site after the spac- 
ing is increased from 2.2 /xm to 5.5 /xm for different ramp 
times. Vo = h X 37 kHz. The inset in (a) shows the excited 
state probability against ramp time for the outermost lattice 
site {ux = 8, n^ = 8). (b) Probability of finding an atom 
in a bound state for the same parameters. The lines denote 
the numerically calculated probability and the points show 
experimental data. 



We perform a more detailed analysis of the minimum 
expansion time by numerically solving Eq. Q. The op- 
tical lattice potential V{x^y^ t) has the form given in Eq. 
(IT]), with our experimental parameters (Vb = hx37 kHz, 
do = 2.2 /im, w = 69 /im, S = 2.5) and the harmonic 
trapping potential is included. We take the initial wave 
function |'0n^,n {t = 0)) to be the ground state of a sin- 
gle well of the sinusoidal potential centered on x = n^^do, 
y = riydo. A 6*^ order finite difference formula is used 
to express the second derivative, and the wave function 
is evolved using the propagator U{t) = ex.-p[—iH{t)At/h] 
for each time-step At. We solve Eq. Q in one dimension 



The probability of an atom remaining in the gr ound 
state, \{^l^n^,ny{t = t^)|V^ground)P, IS showu in Fig. |3(a) 



(we find the probability for the n*^ lattice site by aver 
aging over Uy for each value of Ux). For the outermost 
lattice site {rix = 8, n^y = 8) the ground state fidelity is 
more than 99.99% for ramps times longer than 12 ms. 
The time taken to spontaneously scatter a photon is 7 s 
in our lattice (with Vq = hx 37 kHz), and can be several 
minutes in optical lattices with blue frequency detuning. 
Thus many expansion and addressing operations could 
be carried out with high fidelity before coherence is de- 
stroyed due to spontaneous emission. 



The probability of an atom remaining in a bound state 
of the optical lattice is also numerically calculated and 
we compare this with experimental results. In the ex- 
periment, we expand the lattice as described above for 
different ramp times. To determine the probability of an 
atom in a given site remaining bound, we first sum over 
the columns of the lattice, then calculate the number of 
atoms in site n. The number of atoms in the outermost 
site does not change (within the margin of error) for ex- 
pansion times between 8 and 20 ms, so the number of 
atoms is normalized to the number of atoms in s ite n for 
these ramp times. The results are shown in Fig. |3(b)[ 

The experimental and numerically calculated results 
agree well. The small discrepancy (the ramp time for 
a given probability is slightly longer in the experiment) 
may be due to the fact that the numerical results are 
for a single particle only. In the experiment, however, 
there are on average a few hundred atoms per lattice site, 
and the calculations do not take into account the energy 
due to repulsive interactions. The close agreement with 
the numerically calculated results indicates that the main 
source of atom losses is the speed of expansion rather 
than from a source of instability in the ramp, confirming 
that the ramp of the optical potential is smooth. 

Note that the minimum ramp time has been calculated 
for only one lattice depth and one form of ramp. With 
a deeper lattice, the expansion timescales would be even 
faster. Also, it is possible that the minimum ramp time 
could be reduced further by optimizing the form of the 
ramp. 

In this experiment the initial lattice spacing was 
2.2 /im, and this technique can be easily used to ex- 
pand from smaller lattice sizes if a higher numerical aper- 
ture (NA) lens is used to form the optical lattice — the 
minimum lattice spacing dmin = A^/2NA, where Xi is 
the wavelength of light forming the optical lattice. Note 
that (imin is not equal to the imaging resolution; atoms 



can just be resolved when they are spaced by more than 
dres = 0.6Ai/NA, where A^ is the wavelength of imaging 
light. For A^ < A^ (blue frequency detuning) the lattice 
spacing is smaller than the imaging resolution. 

In our experiments on the accordion lattice we have 
concentrated on proving its usefulness for expansion of 
cold atoms in a periodic potential which is equivalent to 
increasing the magnification of the optical imaging sys- 
tem, however there are other modes of operation of this 
new apparatus: (1) Two, or more, optical lattices of dif- 
ferent spacing can be created at the same time by apply- 
ing the requisite driving frequencies to the AODs. The 
light intensities and resulting optical potentials add to- 
gether forming a superlattice. In future work cold atoms 
could be loaded into an optical lattice of double wells by 
loading atoms into a lattice of period d, expanding the 
spacing of the wells to 2d and gradually switching on an- 
other lattice of spacing d to split each potential well in 
two; (2) the technique can be extended to three dimen- 
sions by applying an additional accordion lattice in the 
axial direction, enabling planes of atoms to be moved out 
of focus of the imaging system, thus facilitating imaging 
of a single plane; (3) the optical lattice can be expanded 
to decrease the tunneling energy without changing the 
lattice depth, allowing investigation of different quantum 
phases. 

The use of an accordion lattice for expanding the pe- 
riod of an optical lattice has been demonstrated for the 
case of absorption imaging, but it could also be utilized 
for single atom fluorescence imaging. The arrangement 
used to generate the optical lattice is extremely flexible 
and allows smooth variation of the optical potential as 
shown here, and demonstrated previously with a rotat- 
ing optical lattice [20] . 
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